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ON THE propagation OF LONG WAVES IN ACOUSTICALLY TREATED, CURVE DUCTS 

by w. RosLdfinsKi 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 441 Jd 
abstract 

A two-dimensional, detailed study is presented on the oehavior of long 
waves in lined, curved ducts. Tne analysis includes a comparison between the 
propagation in curved and straight lined ducts. A parametric study was 
conducted over a range of wall admittance and duct wall separation. The 
complex eigenvalues of the characteristic equation, whicn in the case of a 
curved duct are also tne angular wavenumuers, have been obtained by 
successive approximations. 

INTRODUCTION 

There are both scientific and practical reasons tor understanding the 
mechanics of the propagation of acoustic waves inside acoustically lined 
pipes and oucts. The. noise suppressors found in aircraft engine ducts are a 
typical application. Consequently, many papers and engineering articles have 
been published on the motion of waves in straight conduits with acoustically 
absorbent, or lined, walls. 

There are some instances where extreme curvature is present in an engine 
duct. An example is the S-snaped inlet to tne center engine of a 
three-engine airplane. For some time now, mere has been a theoretical 
effort focused on tiie propagation of waves in lined, curved ducts. Of about 
12 papers published to date on the propagation of waves in circularly curved 
ducts, three treat tne case of soft-wall oucts of rectangular cross section. 
This relative scarcity of literature on this subject may ut explained by the 




severe mathematical difficulties that exist when eigenvalue problems are 
encountered which involve Bessel functions and other higher transcendental 
functions, especially when the functions are complex.. In the three papers, 
selected areas of the propagation process are developed. 

The first paper, by Grigor'yan^, published in 1970, considers the 
numerical solution of the propagation equations at two f i ‘^quenc ies, 2 and 4 
kHz. Three conf igurat ions of acoustical wall treatment are considered: both 

curved walls lined; only the convex wall lined; ana only tne concave wall 
lined. The radii of curvature of the bends are equal to or greater tiian the 
radial distance between the curved walls, so that a = R^/R| < 2, where 
R£ is the radius of tlie concave wall, and is the raoius of the 

convex wall. This is a rather narrow range, considering that in the great 
majority of industrial appl icai ions, a > 2. The analysis considers two 
widttis of tne curved ducts, and two types of sound absoroing materials. A 
number of very interesting conclusions are obtained, among them t..e fact that 
the curvature does not always yield an increase in the attenuation over that 
obtainable with a straigiit duct. 

Myers and Mungur^, in a paper published in l97b, numerically solve the 
partial differential equation of motion, and numerically evaluate the 
propagation ano attenuation of three duct mooes. The zeroth, basic mode is 
not included. In this work, they assume one arbitrary value of wall 
admittance, a range of vibrational frequencies, ana a single duct width 
characterized by a = 2. In Uieir conclusions, the authors state that for 
locally reactive duct walls, the sound field generally attenuates more 
rapidly along the axial direction of a bend than does the field in a straight 


duct. 
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Finally, in 1977, Ko anu Ho"^ puolisi*e<l results of tneir extensive 
parametric stuay of motion in slightly bent oucts. The purpose of their work 
was to develop a method for predicting the attenuation of sound in 
acoustically lined, curved ducts. They proceed by the derivation of 
eigenvalues, which they list in extensive tables. Among their finuings was 
tlie determination that the level of attenuation of sound is rather dependent 
on acoustic resistance and relatively insensitive to tne sharpness of the 
bend, within the limits 1.11 < a < 1.43, only the median arc length affects 
the amount of souno attenuation, which is independent of tne width of the 
duct. 

The present study is intended to extend this author's two-dimensional 

study of the behavior of long waves in circularly bent ducts with hard 
4 

walls. One of the results of that work was determination of wave phase 
velocity, which in the case of very low frequencies is higher than in a 
similar straight duct and is proportional to tne sharpness of the bend. This 
extension of that study consists of modifying the boundary conditions and 
evaluating the effect of the new boundaries on sound phase velocity and on 
the oistribution and level of particle velocity and pressure. These boundary 
conditions assume that the curved walls of the bend will be lined with 
locally reactive, acoustically absorbing material. All other assunptions of 
tne previous work are accepted without cliange, namely that there is a 

sustained, continuous, and steady propagation of long acoustic waves in a 
two-dimensional, infinitely long, circularly bent duct. The closest physical 
example of such an infinitelj’ long bend is a tightly wound coil of which the 
pitch is negligible compared to the radius of curv.uure of the duct. The 
wave length is at least two orders 0 ‘ magnitude larger than the width of the 
duct. The pressure distribution within the duct satisfies the wave equation. 
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ANALYTICAL SOLUTION 
Infinite, Lineu, Curved Duct 

The basic difficulty in solving equations of motion with complex roots 
resides in our inability to split the final equation containing exponential 
and trigonometric functions into real and imaginary parts. Were this 
possible, it still would not guarantee yielding the answers directly, but it 
would help enormously, because a single term would be sufficient to satisfy 
the characteristic equation. With complex expressions, two real numbers must 
be found. So we are forced to rely on the tedious method of successive 
approximations. 

The calculation of pressure and of components of particle velocity starts 
wi tn def i nit ion of the ve locity potentia 1 : 

<t> [A^J^(kr) + B^Y^{kr).] expi(wt - ve) (i) 

V 

which yields 



v 


e 


1 

r 30’ 


and 



It 

3r 


A sympol list is found in the appendix. In the above expressions, the two 
bessel functions J and Y are of complex order, since the angular 
waven umber 


V = V i v^ 

Determinat ion of cotistants in the definition of the velocity potential 

is usually obtained from the boundary condition, by assuming tnat at e = 0, 
V. = V . However, in this paper only the lowest order inod'e will be 
considered. The attenuation of sound, in decibels (dB), is defined by adB = 
20 • Since p ^ exp(-v^e), adli = 20 (-v^g) log e. The 

angular wave phase velocity, o = u/v^, and tne phase velocity along tne 


5 


median radius, s = = c kRj^(a + l)/2vp are easily calculated. A 

general solution of the two-dimensional problem for both rigid and soft-wall 

»v 

cases has been formulated by Grigor'yan and others. Tne equation for 
determination of the angular wavenumber for the case of identical lining 
material on both curved walls (the original equation makes it possible to 
stipulate different materials) of a circularly bent duct, is 


[j'(kR^)Y^(kR2) - J^(kR^)Y‘(kRj)j + ir,[j‘(kR^)Y^(kR^) 

- J^(kR2)YjkRi) - J^{kRi)Y‘(kR.O ^ j’(kR^)Y^(kR^)] 

+ /[J^(kRpY^(kR^) - J^(kR2)Y^(kRp = 0 (2) 

where n = -pc v^./p (i.e., a dimensionless wall aoiiiittance v,./p). In 
general, n - t + io is complex, and t is the conductance of the walls, 
and 0 tile susceptance of the walls. It is clear that using r, = 0 in 
Eq. (2) leaves only tne expression of the cross products of the derivatives 
of the two ISessel functions, which corresponds to tne case of the hard-walleu 
bend. 

To solve Eq. (2) for the case of ‘long waves propagating in a soft-walled 
bend, the two Bessel functions are expanded in series, and in view of the 
smallness of the arguments, only the first terms of the series are retained. 
Tne second and following terms of the expansions are smaller. than tlie first 
term at least by the order of (kR^ . Rearranging and eliminating small 
terms of higher oraer yields 


2 

T 


2iro (< 


it) 


a + 1 



2 

V 


a( kRj^) 



For T = o = 0, the conditions for stiff walls, Eq. (3) yields a'* = a~'' 
or a^'’ = 1, and v = inn/ In a, pure imaginary roots derived before in 


( 3 ) 
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Ref. 4.. Eq. (3) also allows verification that for finite values of t and 
a, there are no pure imaginary or simple real roots. Hence, v must be 
complex. Roots, which are complex angular wavenumbers, have been obtained by 
successive approximations. It helped to know that both tne real and the 
imaginary parts of u must be small oecause tlie real root for the case of 
the hard walls is small. In applying Eq. (3), calculations have been made 
for the following set of parameters: k = 0.1; = 0.2; a = 2 and 4; 

T = 0.1, 0.2, ..., O.b; and o = 0.2, 0.3, 0.6, and 0.8. The calculated 
complex roots of Eq. (3), which are also complex propagation constants, are 
listea in Table I and are shown in Fig. 1. It will be noted ttiat quantities 
related to attenuation are one oraer of magnitude smaller than terms related 
to propagation. Attenuation increases with increasing t, but attenuation is 
less pronounced witli a wider duct. For small values of conductance t, 
susceptance o does not greatly influence attenuation; with higher t, 
however, large changes occur with changing o. The percentage change in 
attenuation at constant t and varying o is nearly twice as great for the 
lower range of t as for the largest t ■ studied.. On the same graph, the 
angular wavenumber for tne rigid, curved wall case (i.e., for t = o = 0, 
and a = 2) is also shown. Presence of acoustical lining on curved walls 
markedly increases ttie values of the propagation constants 

Infinite, Lined, Straigiit Duct 

The propagation of acoustic waves in a straight, lined duct can be 
described in tenns of the acoustic pressure using the expression given by 
Rice^: 

p = cos ~ c 

In Eq. (4), the'g's are complex roots. Also, ^ = a + is is a complex wave 
number with its propagation and attenuations terms. Now, f = - (2gykH)2 
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as derived by Rice, H/2 is the half-width of the duct, and x and y are 
the axial and transverse coordinates, with y = 0 at duct centerline. 

The phase velocity in a straight duct will be simply x = c/B. The 
transverse particle velocity may be expressed by u^ = nP/pc or 
= i{ 3p/sy)ypu). Equating these two expressions at y = H/2 yields 

“"'9> - ■ (f ?)(i) 

wiiere o = t+ io, as in the case of the curved uuct. Solution of Eq. (5) 
yields the real and the imaginary components of the roots. Following this, a 
calculation of complex wavenumber is straightforward. Also, the expression 
for the axial component of particle velocity is easily obtainaoie by 
derivation of the expression for pressure with respect to x, the axial duct 
coordinate [u^ « i(<ip/3x)/phi]. 

To be able to compare the propagation in.tne straight duct with the 
propagation in the curved duct, not only the same set of values for r and 
a are used in numerical calculations as in the case of the curved duct, out 
also the widths of the straight ducts are selected to match the two widths 
used ih analysis of the curved ducts. The complex wavenumber ? for the 
straight, lined duct (i.e., the attenuaton and the propagation terms) are 
given in Table II and shown in Fig. 2. As in the case of a curved, lined 
duct, the wider duct has a less pronounced attenuation. The general 
influence of t and o on the complex wavenumber is similar to that in the 
curved duct. 

RESULTS AND DISCUSSION 
Phase Velocities 

Using ttie propagation constants for both curved and straight ducts, 
corresponding phase velocities were calculated and ratios s/x established. 
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They are shown on Fiy. 3 as a function of t. The ratios are greater than 1 
for all parameters considered. They are higher for a wider duct (e.g., for 
a = 4) than for one half as wide (for a = 2). Generally, the phase 
velocities and their ratios increase with t, but not strongly. The 
corresponding values of phase-velocity ratios for a hard-walled oucts are 
also shown. At these low frequencies, the phase-velocity ratios in lined 
ducts are lower than in untreated ducts. 

Attenuation of Sound 

Attenuation of sound in lined ducts, both curveo and straight, has been 
calculated for all parameters shown in Tables I and II. Some of the results 
for curved ducts are shown in Fig. 4, in whicn pressure oecay in a bend of 
180° is plotted against parameters of the lining, both t ana o. 

There is a very significant increase in pressure loss with the increase 
of conductance t. On the other hand, for a given level of conductance, 
higher values of susceptance o somewnat reduce the pressure loss. Both 
relations are clearly nonlinear. Also, it is ciaar that in a wider bend the 
losses are somewhat less than in a narrow one. It will be noticeo that on a 
plot of pressure ratio versus t, the two curves can be extrapolated down to 
T = 0 (i.e., to the case of bends with haro walls). The extrapolation 
reaches a pressure ratio of unity, ar it should. 

Results of calculations of sound attenuation per unit length, in dB/m, 
ana comparisons of attenuation between curved ana straight ducts, are shown 
in Fig. 5. Attenuation has been evaluated witn variable parameters, t and 
o, of the acoustical lining in ducts of two wiaths. Sound attenuation 
rapidly increases with increasing t, but it is. less pronounced in a wider 
duct. In general, attenuation is slightly higher in a straight duct (by 2 to 
2 Percent), but in this range of very Jong waves, it must be considered 
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significant. All four linear functions shown in Fig. 5 extrapolate to zero 
attenuation at t =0. The dependence of attenuation on o is nonlinear, 
and attenuation decreases with increasing susceptance. For narrow ducts, the 
difference between the attenuation in a straight duct and in a curved duct 
remains fairly constant with increasing o. For wide ducts, this difference 
in attenuation uecreases with increasing o. 

Radial Distribution of Wave-Motion Parameters 
The radial distribution of the particle velocity components and of the 
acoustic pressure has been analyzed for all cases taken into consideration 
previously. Samples of typical profiles in a curved, lined duct are shown in 
Fig. 6. The particle velocities are nondimensional ized by using v^. 
Basically, very slight radial variations have been detected when t and o 
were given their full range of values. Also, when the wave was moving oown 
the curved duct, the pressure and axial-velocity pr^'^iles remained unchanged 
except for a gradual decrease in their amplitudes. The profile of the, radial 
component of the particle velocity is different in, that its slope changes 
with distance. The radial velocity distribution exhibits a zero point near 
the convex wall of the bend. This point isriot much displaced toward the 
center of the duct when a duct three times as wide is used. In Fig. 7, these 
profiles of particle velocity components are compared with distrioutions 
calculated for the case of unlined bends of the same geometry and for 
similar, very low frequencies of acoustic waves. A striking di'ference 
exists between the two radial particle velocity distributions. This vvas to 
be expected, because a lined duct has a finite value of wall impedance wnich 
does not require tnat the radial particle velocities vanish at the wall. The 
relative values of the radial velocities show; on the graph have no 
particular meaning, because the boundary conditions at e = 0 for the rigid 
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■i and the lined bends are different. The distribution of the tangential 

: velocity component is not affected by the presence of a lining, except for a 

generally lower amplitude at all radii. 

The radial distribution of the acoustic pressure is markedly changed by 
! the presence of a lining. The data are nonaimensinal ized by using p^, the 

[ ■ acoustic pressure corresponding to reference particle velocity v^. within 

; a hard-walled bend, the pressure is cnaracteristica 1 ly higher at the outside 

‘ 7 

i wall, as was first documented by Cunmings . In a lined, curved duct, it is 

i practically independent of the radial postion. 

i 

' CONCLUDING REMARKS 

' Propagation of long waves in circularly bent ducts lined with 

; sound-absorbing material has been analyzed for a range of acoustic 

■ parameters. A set of eigenvalues, the complex angular wavenumberS; has been 

‘ obtained by a method of successive approximations. The results have been 

compared with the propagation of low-frequency sound in a straight, lined 
I : duct of the sane width and lined with the same material as the bent duct. 

I . . • ■ 

I ‘ Furthennore, a comparison witn the propagation parameters in a hard-walled, 

; bent duct has been presented as well. 

At these very low frequencies, and in the range of parameters evaluated, 

, the sound attenuation in a curved, lined duct is generally less pronounced 

; than in a straight duct of ttie same width and same wall impedance, 

r Results of this study and of work of others indicate that there is need 

-r for furtiier analysis of the propagation of sound in curved, lined ducts. The 

; overall picture of such propagation is not yet clear. 
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APPENDIX - LIST OF . SYMBOLS 

a constant 
radii ratio = R 2 /R 1 
a constant 
velocity of sound 

complex eigenvalue for the straight duct 

depth or width of a rectangular duct 

Bessel function of the first kind, of order 

wavenumber 

acoustic pressure 

reference acoustic pressure 

acoustic pressure at angular location e 

curved-duct convex wall radius = 0.2 m 

curved-duct concave wall radius 

radial coordinate 

distance along curved-duct centerline 
wave phase velocity in the curved duct 
time coordinate 
reference particle velocity 
tangential component of particle velocity 
radial component of particle velocity 
axial component of particle velocity 
transverse particle velocity 
axial coordinate ' 

wave phase velocity in a straight duct 
Bessel function of the second kind, of order 
transverse coordinate 






i o attenuation coefficient in a straight duct 

I 

i 8 propagation coefficient in a straight duct 

; t complex wavenumber, t = a + iB 

n wall admittance 

e polar coordinate 

0 angular phase velocity 

V complex angular wavenumber, v = + iv2 

propagation coefficient 

v2 attenuation coefficient 

p density 

a susceptance 

T conductance 

0 ' velocity potential 

u angular frequency = kc 
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TABLE I - COMPLEX ANGULAR WAVENUMBERS FOR CURVED LINED DUCTS FOR 


kRi o 0.02 AND FOR A RANGE OF PARAMETERS OF THE LINING 


T 


a = 

2 

a = 

4 

Propagation 

''1 

Attenuation 

''2 

Propagation 

''1 

Attenuation 

"2 

0.1 

0. 

2 

0.13S46 

-0.03201 



.2 



. 144bO 

-.05997 



.3 



.ibb?; 

-.08351 

0.14225 

-0.07680 

.4 



.Ib73b 

-.10364 



- .S 

■ ' 


- .17873 

-.12131 



.1 

• 

5 

.20931 

-.02077 

.19199 

-.01925 

.2 



.21224 

-.04097 

.19462 

-.03797 

.3 



.21o72 

-.06019 

,. 19864 

-.05580 

.4 



.22234 

-.07823 

.2036/ 

-.07256 

.b 



.22873 

-.09505 

.20938 

-.08823 

.1 

• 

8 

.26415 

-.Olbbi 



.2 



.26565 

-.03283 



.3 



.26805 

-.04860 

.24649 

-.04564 

.4 



.27125 

-.06480 




.5 

ij 

) 

.27511 

-.07925 
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TABLE II. - complex WAVENUMBERS FOR STRAIGHT LINED DUCTS FOR 


k = 0.1 AND FOR A RANGE OF PARAMETERS OF THE Li.wNG 


■ 

1 

H72 

= 0.1 

H/2 = 

0.3 

Propagation 

ktj 

Atteruialion 

ko 

Propagation 

kB 

Attenuation 

ko 

0.3 

0. 

2 

0.52195 

0.28687 

0. 30218 

0.18561 

.1 


5 


.07175 

. .39872 

• .05574 

.2 



.71197 

.13971 

.40343 

.08171 

.3 



. 72886 

.20486 

.41296 

.11986 

.4 



.74772 

. 26620 

.42498 

.15527 

.5 

1 

f 

.76908 

.32183 

.43856 

.18817 

.3 


S 

.90287 

.16522 

.51268 

.08939 
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Fig. 1. - Real and imaginary parts of complex angular wavenumbers fora 
curved duct 




PROPAGATION COEFFICIENT, P 

g. 2. - Real and imaginary parts of complex wavenumbers for a 
straight duct. 
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(al PARTICLE VELOCITY COMPONENTS AT 0 = u FOR TWO 
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(b) PARTICLE VELOCITY DISTRIBUTION CHANGES WITH DISTANCE. 
FOR T = 0. 1, 0 = 0. 5. AND FOR TWO VALUES OF a. 

Fig. 6. - Particle velocity components in a lined curved duct 
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